Chinay ww 会 1 在 甘 9 车 || 
ChinaA IV r1 F BH ti 


大 口径 射电 望远镜 轮 轨 滚动 接触 有 限 元 分 析 
师 洪 强 !2 ， 许 谦 234， 王 娜 234， 王 兆 军 1 
(1. 新 疆 大 学 物理 科学 与 技术 学 院 ， 乌 鲁 木 齐 ，830046; 
2. 中 国 科学 院 新 疆 天 文 台 ， 乌 鲁 木 齐 ，830011; 
3. 中 国 科 学 院 射 电 天 文 重点 实验 室 ， 乌 鲁 木 齐 ，830011; 
4. 新 疆 射 电 天 体 物理 重点 实验 室 ， 乌 鲁 木 齐 ，830011) 
摘要 : 通过 显 式 有 限 元 法 建立 轮 轨 滚动 接触 模型 ， 用 于 分 析 110 米 口径 射电 望远镜 滚轮 
经 过 焊 颖 时 的 应 力 变 化 。 建 模 时 不 仅 考虑 了 焊 颖 处 材料 差异 ,同时 还 考虑 轨 道 焊 颖 处 表面 不 
平 度 和 摩擦 系数 的 影响 , 分 析 了 滚轮 经 过 轨道 焊 颖 时 轮 轨 间作 用 力 的 瞬 态 响应 采用 层级 细 
分 技术 对 轮 轨 模型 进行 前 处 理 ， 经 过 静 力学 分 析 结 果 表 明 ， 模 型 求解 精度 较 高 。 将 该 前 处 理 
模型 用 于 瞬 态 分 析 发 现 , 焊 颖 处 表面 不 平 度 对 轮 轨 间 作用 力 影响 较 大 , 使 轮 轨 间 Mises 应 力 
增加 了 约 20%， 可 能 会 对 轨道 长 期 可 靠 服 役 带 来 一 定 的 隐患 ， 因 此 需要 定期 对 轨道 的 状态 
进行 检测 与 维护 ， 保 障 轨道 的 长 期 可 靠 服 役 特性 。 
关键 词 : 射电 望远镜 ; 方位 轨道 ; 滚动 接触 ; 瞬 态 响应 ; 有 限 元 法 
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中 小 口径 精度 要 求 不 高 的 射电 望远镜 普遍 采用 拼接 轨道 技术 , 拼接 轨道 在 滚轮 经 过 接 缝 
时 易 出 现 跳动 以 及 应 力 集中 ， 从 而 影响 天 线 指向 精度 以 及 轨道 使 用 寿命 山 。 目 前 大 口径 、 高 
精度 射电 望远镜 轨道 逐渐 采用 整体 焊接 技术 , 如 美国 GBTPI、 上 海天 马 65mB] 射 电 望 远 镜 等 。 
正在 建设 的 新 疆 110 米 口 径 全 向 可 动 射 电 望 远 镜 (QTT) 内 ， 方 位 轨道 将 采用 整体 焊接 技术 ， 


- 由 于 采用 窗 档 坡 口 焊接 技术 ， 因 此 堆 焊 层 需要 用 到 3 种 不 同 的 焊接 材料 四。 焊接 过 程 中 的 高 


温 会 引起 材料 内 部 结构 性 能 变化 , 造成 焊 颖 位 置 硬度 和 强度 与 母 材 的 差异 。 由 于 轮 轨 接触 力 
分 布 有 大 载荷 、 小 集中 的 特点 ， 天线 深 轮 在 经 过 焊 颖 位 置 时 ， 焊 颖 与 母 材 的 差异 可 能 会 引起 
轮 轨 接 触 力 的 突变 。 应 力 突变 是 引起 疲劳 磨损 的 主要 原因 之 一 , 因此 弄 清 深 轮 经 过 焊 颖 时 的 
应 力 变化 尤为 重要 。 

对 于 大 型 射电 望远镜 轨道 而 言 ,传统 试 验 研 究 成 本 高 昂 , 因此 可 采用 仿真 分 析 的 方法 进 
了 研究 中 。 任 就 等 人 中 针对 大 型 反射 面 天 线 轮 轨 接 触 问题 , 在 不 同 仿真 软件 中 进行 对 比分 析 ， 
得 轮 轨 接触 应 力 ， 但 未 涉及 瞬 态 分 析 。Juneja GIB 等 人 通过 三 维 有 限 元 分 析 ， 用 移动 载荷 
法 模拟 轮子 滚动 , 分 析 GBT 和 斜 接 颖 轨道 的 力学 行为 以 及 耐 磨 板 、 基 板 接 触 界面 的 微 动 磨损 ， 
分 析 结 果 为 轨道 重新 设计 与 改造 提供 参考 。 移动 载 苟 法 不 能 还 原 真 实 的 轮 轨 接 触 状态 , 且 材 
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料 差 异 带 来 的 影响 也 很 难 考 虑 进去 。 因 此 需要 采用 完整 的 有 限 元 模型 才能 更 加 真实 、 准 确 的 
反映 轮 轨 接触 的 瞬 态 响应 。 

本 研究 采用 有 限 元 法 将 天 线 轨道 焊 缝 真实 焊接 工艺 考虑 在 内 , 建立 QTT 轮 轨 接触 模型 ， 
使 用 层级 细 分 技术 对 模型 进行 相应 处 理 。 在 此 基础 上 采用 “ 面 - 面 ”赫兹 接触 算法 求解 接触 
而 合 , 进而 通过 显 式 时 间 积 分 法 , 获得 轮 轨 间 作用 力 在 不 同 摩擦 系数 和 轨道 面 形 貌 下 的 瞬 态 
响应 。 研 究 结果 为 天 线 轨 道 的 疲劳 寿命 分 析 和 天 线 维护 提供 参考 。 


1 轮 轨 接触 应 力 计算 


本 研究 建立 的 大 型 射电 望远镜 轮 轨 接触 模型 , 由 于 接触 斑 的 几何 尺寸 远 小 于 轮 轨 几何 特 
征 尺 寸 , 且 不 考虑 天 线 滚轮 的 侧 滑 和 里 滑 现 象 。 因 此 轮 轨 间作 用 力 可 以 使 用 赫兹 接触 理论 进 
行 计算 匀 ， 通 过 理论 计算 与 仿真 结果 进行 对 比 ， 验 证 模型 的 准确 性 ， 为 瞬 态 分 析 黄 定 基础 。 

根据 Hertz 接触 理论 有 关 弹 性 体 接触 斑 的 几何 尺寸 关系 ， 轮 轨 之 间 的 接触 可 看 作 是 两 个 
任意 曲面 弹性 体 的 接触 ， 它 的 接触 斑 呈 椭圆 形 ， 其 长 轴 半 径 a， 和 短 轴 半 径 5 分别 为 : 


men 3zN (k, +k) v 3z N (k, *- k,) (D 
4(A-* B) 4(A-* B) 
N 是 轮 轨 间 的 法 向 力 ; m 和 与 椭圆 积分 相关 ， 一 般 表 示 为 : 
Ry (2) 
A+B 
m. n. 0 的 关系 可 以 通过 查 表 获 得 。h 和 态 为 材料 参数 ， 其 值 表示 为 : 
1—w _1-v; 


(ok 
7 E, zE, 


其 中 ，B1、E2 和 vj、w 分 别 为 滚轮 和 轨道 材料 的 弹性 模 量 和 泊 松 比 。 轮 轨 接 触 曲面 初始 
间隙 函数 的 常数 4 和 B 可 由 式 (4) 确定 。 


acre) 2E (4) 
2 Ry Ry R4 Ry 2 Ri R R4 Ry 


k, = (3) 


式 中 Rjy、R12、R21、Rz2 分 别 为 两 个 接触 体 的 主 曲 率 半 径 ， 求 得 接触 斑 长 、 短 半 轴 a、b 
之 后 ， 则 接触 斑 的 内 应 力 分 布 即 可 表示 为 : 


3 N xy y ? 
PAS ah i B B en 
则 由 式 〈5) 可 知 ， 最 大 接触 压 应 力 g 出 现在 接触 斑 中 心 处 ， 其 值 为 : 
3N 
一 (6) 
do 2zab 


由 于 接触 位 置 多 为 压 应 力 ， 轨 道 一 般 不 会 因此 破坏 ， 而 如 果 剪 应 力 超 过 容许 值 ， 会 使 得 
轨道 内 部 发 生 塑性 流动 , 更 容易 引起 钢轨 的 破坏 。 最 大 剪 应 力 发 生 在 轮 轨 接触 面 以 下 的 某 一 
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深度 ， 其 值 约 为 : 


PE 
2T ~ 0.63q, = 0.63m, (7) 


2 
R; 


接触 面 上 的 最 大 前 应 力 为 : 
27, = nq, (8) 
式 中 mo、n0 与 两 接触 体 的 主 曲 率 半径 有 关 ， 均 可 由 查 表 得 到 09。 通 过 以 上 公式 即 可 计 
算出 天 线 轮 轨 接 触 力 的 大 小 和 分 布 。 


2 轮 轨 接触 模型 构建 


2.1 模型 建立 

仿真 模型 基于 QTT 建立 ， 天 线 座 架 共有 4 组 滚轮 ， 每 组 8 个 。 取 一 个 滚轮 和 两 段 轨道 
组 成 分 析 模 型 ， 如 图 1 所 示 。 轨 道 建 模 时 考虑 了 焊接 过 渡 层 、 填 充 层 及 表面 焊接 层 。 滚 轮轴 
采用 3 维 线性 有 限 元 应 变 梁 单元 表示 , 并 使 用 多 点 约束 单元 将 轮轴 与 天 线 滚 轮 进行 耦合 ,不 
仅 可 以 将 整个 天 线 载荷 均匀 施加 在 滚轮 上 , 也 可 以 实现 滚轮 的 滚动 , 有 效 地 简化 了 模型 又 不 
失 准 确 性 。 

有 限 元 模型 的 计算 准确 性 高 度 依赖 于 网 格 划分 精度 , 当 接 触 斑 区 域 网 格 尺寸 为 接触 斑 短 
半 轴 宽度 的 1/20 时 ， 可 以 得 到 精确 的 轮 轨 接触 解 ， 当 网 格 尺 寸 为 110 时 ， 可 以 满足 大 多 数 
工程 问题 中 的 精度 要 求 t0。 采 用 层级 网 格 细 分 技术 ,对 焊 颖 附近 区 域 的 轨道 及 滚轮 进行 网 格 
细 化 ， 最 小 网 格 尺 寸 为 1.1 mm。 如 图 1 所 示 ， 模 型 采用 8 节点 6 面体 单元 进行 离散 ， 焊 颖 
位 置 采 用 共 节 点 方式 ， 以 减少 接触 设置 ， 提 高 计算 效率 。 为 了 保证 网 格 质量 ， 所 有 网 格 的 长 
宽 比 小 于 3.5， 最 终 整 个 模型 网 格 数量 约 79 万 ， 节 点 数 82 万 ， 其 中 焊 颖 附近 网 格 数量 约 占 
总 体 网 格 数量 的 91%. 


图 1 模型 及 网 格 


Fig.1 The model and grid 


轨道 材料 为 高 质 合金 钢 42CrMo， 而 焊 颖 位 置 表面 需要 较 大 硬度 以 保证 其 耐 磨 性 和 使 用 


寿命 ,因此 焊接 坡 口 需要 一 层 较 “ 软 ” 材 料 作为 过 渡 ， 将 轨道 材料 与 表面 堆 焊 层 进行 有 效 衔 
接 。 过 渡 层 厚度 约 4mm， 由 于 目前 有 限 元 分 析 无 法 实现 对 材料 硬度 的 表征 ， 因 此 考虑 通过 
采用 不 同 的 弹性 模 量 来 近似 表征 不 同 焊 接 工艺 的 熔 甫 金属 硬度 。 文 献 [12] 通 过 实验 分 析 指 出 
越 靠 近 焊 缝 位置 ， 其 材料 弹性 模 量 越 小 ， 焊 颖 中心 位 置 弹 性 模 量 最 小 ， 约 为 母 材 的 9096. 
因此 取 过 渡 层 及 填充 层 的 弹性 模 量 约 为 轨道 的 900%， 表 面 堆 焊 层 约 为 轨道 的 95%。 具 体 材 
料 参 数 设置 如 表 1 所 示 。 


表 1 材料 参数 


Tablel Material parameter 


Parameter information Elasticity modulus (GPa) Poisson's ratio Material density (kg/m?) 
Rail and wheel 213 
Filling weld 190 
0.3 7850 
Backing weld 190 
Surface weld 200 


2.) 载荷 与 工 况 设置 
考虑 到 天 线 的 重量 以 及 轨道 的 合理 承载 , 将 单 轮 载荷 设 为 200 吨 , 保证 实际 工 况 安全 性 。 
载荷 均 布施 加 于 轮轴 上 , 考虑 轨道 自重 影响 其 底面 设置 为 全 约束 ,滚轮 约束 横向 自由 度 以 避 
免 横 向 滑 移 。 
工 况 一 : 当 滚 动 速 度 为 20 mm/s， 摩 擦 系数 分 别 为 0.3、0.4、0.5 时 ， 分 析 当 天 线 滚轮 经 
过 焊 颖 时 ， 轮 轨 间 作用 力 的 瞬 态 响应 。 

工 况 二 : 当 滚 动 速度 为 20 mms, 摩擦 系 数 为 0.3， 分 析 焊 缝 表面 不 平 度 分 别 为 0.1 mm. 
0.2 mm、0.3 mm 时 ， 对 轮 轨 间作 用 力 的 影响 。 

瞬 态 分 析 时 ， 设 定 滚轮 与 轨道 接触 起 始 位 置 与 焊 颖 轨道 交界 线 之 间 为 稳 态 过 渡 区 ， 以 保 
证 滚轮 进入 焊 缝 区 时 初始 激 扰 的 能 量 被 消耗 031。 通 过 试 算 ， 稳 态 过 渡 区 取 值 为 70 mm. 
2.3 时 间 步 长 确定 

时 间 步 长 与 数值 分 析 误 差分 布 有 关 , 合理 的 时 间 步 长 可 以 保证 数值 分 析 稳 定性 和 计算 精 
度 。 时 间 步 长 的 选取 不 能 一 次 定论 , 应 该 对 比 三 个 相距 一 定时 间 间 隔 步 长 下 的 计算 结果 ， 若 
相对 偏差 较 小 ， 即 可 取 中 间 那 个 步 长 的 计算 结果 作为 问题 的 合理 解 多 。 取 三 个 时 间 步 长 分 
别 计算 2 s 内 深 轮 经 过 相同 位 移 时 ， 绘 制 轨道 上 同一 位 置 的 应 力 时 变 曲线 ， 如 图 2 所 示 。 从 
图 中 可 以 看 出 ， 不 同时 间 步 长 对 分 析 结 果 影 响 较 小 ， 因 此 选取 时 间 步 长 为 0.025 s 进行 瞬 态 
分 析 。 


二 


Time step 
一 “一 0.01s 
—e— 0.025s 

m 300 zs ss 

z 

$ 200 

9 

三 100 

0.0 0.5 1.0 L5 2.0 


Time / s 
2 时 间 步 长 对 结果 的 影响 


Fig.2 The effect oftime step on the results 


3 轮 轨 接 触 有 限 元 分 析 


3.1 静 力 学 分 析 
基于 赫兹 接触 理论 ， 考 虑 天 线 深 轮 处 于 焊 颖 处 的 静态 工 况 ， 设 置 静 摩擦 系数 为 0.5。 通 
过 有 限 元 数值 分 析 ， 静 态 轮 轨 接 触 总 变形 量 为 0.147 mm， 轨 道 最 大 变形 量 为 0.069 mm, H 
道 受 力 云图 如 图 3 所 示 。 由 图 3(Ca) 可 知 , 天 线 滚轮 与 轨道 间 的 Mises 应 力 最 大 值 为 363 MPa, 
与 理论 值 371 MPa 误差 仅 为 2.2%， 最 大 Mises 应 力 出 现在 轨道 表面 下 约 4mm 处 ; 由 图 3 
(b) 可 以 看 到 轮 轨 接触 广 呈 椭圆 形 ， 接 触 广 宽度 为 13 mm 理论 值 13.2 mm)， 与 赫兹 接触 
理论 基本 吻合 ， 最 大 接触 压力 为 486 MPa CHERIE 489 MPa)， 远 小 于 轮 轨 材料 42CrMo 的 
屈服 应 力 930 MPa， 验 证 天 线 静 态 工 况 下 的 安全 性 ; 由 图 3 (c) 可 以 看 到 最 大 剪 应 力 为 117 
MPa《〈 理 论 值 115 MPa)， 但 远 小 于 轮 轨 材料 的 最 大 前 应力， 且 剪 应 力 场 出 现在 轨道 表面 下 ， 
可 反映 真实 轮 轨 接触 应 力 场 。 静 力学 分 析 验 证 了 模型 的 准确 性 ， 可 用 于 瞬 态 分 析 。 
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J3 (a) Mises 7] zd; (bo 接触 压力 云图 ; (ce) 最 大 剪 应 力 云图 


Fig.3 (a) The pattern for Mises stress; (b) The pattern for contact pressure; 


(c) The pattern for maximum shear stress 


3.2 WEG 


由 于 整个 轮 轨 模型 网 格 数量 较 大 ， 瞬 态 分 析 计 算 量 庞大 ， 对 计算 机 性 能 要 


117 


求 极 高 。 基 于 


模型 中 轮 轨 接触 特点 ， 采 用 减 缩 Z 向 模型 比例 以 减 小 模型 计算 量 来 进行 瞬 态 分 析 。 分 别 选 
W ZI 0.5. 0.2. 0.1 的 比例 模型 ， 其 静 力 学 计算 结果 与 完整 模型 进行 对 比 ， 结 果 如 图 4 所 


示 。 从 图 中 可 以 看 出 ， 随 着 Z 向 模型 比例 的 减 小 ， 减 缩 模型 结果 与 原 比例 模 


型 结果 间 的 误 


差 逐渐 增 大 ， 为 了 兼顾 计算 精度 与 计算 量 ， 最 终 选 取 Z 向 厚度 为 0.2 的 模型 进行 
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图 42Z 向 不 同比 例 模型 计算 结果 


Fig.4 Calculation results of different scale models in Z direction 


3.2.1 摩擦 系数 影响 分 析 


瞬 态 分 析 。 


不 同 摩 擦 系数 会 引起 滚轮 与 轨道 间作 用 力 的 变化 。 摩 擦 系数 分 别 取 0.3、0.4、0.5 与 不 
考虑 焊 缝 处 材料 影响 的 分 析 结 果 进 行 对 比 ， 其 轮 轨 间作 用 力 的 时 变 曲线 如 图 5 所 示 。 轮 轨 间 


Mises 应 力 出 现 两 个 快速 变化 的 峰值 ， 其 原因 是 爆 颖 过 渡 层 与 轨道 存在 变形 差 ， 在 深 轮 接触 


到 过 渡 层 时 ,在 轨道 上 出 现 应 力 集中 ,因此 过 渡 层 承受 载荷 快速 减 小 ,由 于 过 渡 


层 只 有 4 mm; 
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滚轮 将 迅速 转 过 并 抵达 高 硬度 焊接 区 , 应 力 集中 又 出 现在 表面 焊接 层 。 滚 轮 从 轨道 段 进 入 焊 


颖 过 程 中 过 渡 层 承受 载荷 最 小 ， 可 保障 焊 缝 承载 区 域 的 焊接 成 功率 。 

最 大 变化 量 约 22 MPa， 该 应 力 的 突变 在 轨道 长 期 服役 过 程 中 容易 造成 焊 缝 区域 的 疲劳 。 
不 同 摩擦 系数 对 轮 轨 间作 用 力 影响 较 小 ， 如 摩擦 系数 为 0.3 与 摩擦 系数 为 0.5 之 间 的 

Mises 应 力 差 约 为 2.5 MPa， 因 此 摩擦 系数 对 轮 轨 间 作用 力 的 影响 几乎 可 以 忽略 不 计 。 摩 擦 


系数 主要 影响 轮 轨 之 间 的 摩擦 力 。 从 图 5 O) 可 以 看 上 


整个 仿真 中 Mises 应 力 


上 ， 随 着 摩 探 系 数 的 增 大 ， 轮 轨 间 静 、 


动 摩擦 应 力 都 随 之 增 大 ， 且 材料 越 “ 软 ”摩擦 应 力 越 大 。 考虑 焊 缝 影响, 深 轮 在 经 过 焊 颖 时 ， 
最 大 摩擦 应 力 增加 了 71%， 也 增加 了 焊 缝 位置 长 期 服役 过 程 出 现 疲 劳 磨 损 的 可 能 性 。 大 型 
轮 轨 式 天 线 由 于 滚轮 负载 大 、 速 度 低 ， 一 般 不 会 出 现 滑动 摩擦 和 滚轮 “打滑 ?现象 ， 因 此 天 
线 的 日 常 维护 可 以 通过 涂抹 润滑 油 来 降低 轮 轨 间 摩擦 系数 ， 从 而 降低 摩 探 应 力 。 
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到 5 不 同 摩擦 系数 的 影响 。(a) Mises 应 力 时 变 


(b) 


2X; (b) 摩擦 应 力 时 变 曲线 


Fig.5 The influence of different friction factors. (a) The time-varying curve of Mises stress: 


(b) The time-varying curve of frictional stress 


3.22 轨道 表面 不 平 度 的 影响 


QTT 建成 后 整个 轨道 表面 不 平 度 峰 峰值 优 ] 


F0. mm， 而 轨道 在 


程 中 ， 不 可 避免 会 造成 轨道 表面 的 不 平 度 。 由 于 


不 平 度 的 随机 性 和 复 


造 及 现场 焊接 装配 过 
杂 性 ， 使 用 文献 [15] 的 


方法 ， 用 正弦 波 曲 线 简 化 焊 颖 处 表面 不 平 度 。 以 焊 颖 中 心 位 置 为 原点 ， 焊 颖 宽度 为 60mm, 


依照 图 6 D) 所 示 的 正弦 波 曲线 改变 焊 缝 位 置 处 的 节点 匀 


标 以 实现 有 限 元 模型 的 修改 。 


在 不 同 的 表面 不 平 度 下 ， 轮 轨 间 作用 力 的 时 变 曲线 如 图 6 (bo 所 示 。 从 图 中 可 以 看 出 ， 
滚轮 在 进入 焊 缝 时 ，Mises 应 力 出 现 快 速 的 减 小 又 增 大 的 现象 ， 是 因为 通过 坐标 改变 轨道 表 
面 不 平 度 时 ， 曲 面 与 平面 过 渡 处 相 邻 两 个 节点 出 现 了 高 度 差 。 表 面 不 平 度 为 0.1 mm 时 ， 该 
高 度 差 约 为 0.005 mm， 但 造成 应 力 的 剧烈 变化 。 侧 面 反映 了 重 载 天 线 轮 轨 接 触 时 ， 局 部 微 
小 缺陷 (表面 夹杂 、 气 泡 等 ) 更 容易 对 天 线 滚 轮 及 轨道 表面 造成 严重 影响 。 随 着 轨道 表面 不 


4 


平 度 的 增加 ，Mises 应 力也 随 之 增 大 。 当 轨道 表面 不 


F 度 为 0.3 mm 时 ， 最 大 Mises 应 力 为 


442 MPa， 相 比 于 平顺 轨道 其 应 力 增 加 了 约 20%。 可 以 发 现 焊 颖 不平 度 对 天 线 的 影响 较 大 ， 
护 ， 保 证 天 线 轨道 表面 清洁 ， 及 时 处 理 磨耗 


天 线 运行 时 ， 需 要 时 常 对 天 线 轨 道 进行 检测 和 维 
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或 装配 带 来 的 轨道 表面 不 平 度 ， 减 小 对 天 线 轮 轨 的 损伤 。 
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K[ 6 a) 轨道 焊 颖 不 平 度 ;(b) 不 同 轨道 表面 不 平 度 下 的 应 力 时 变 曲线 


Fig.6 (a) The rail welded joint flatness; (b) Time-varying stress curves in different track surface flatness 
分 析 可 知 , 摩擦 系数 对 轮 轨 间作 用 力 影响 较 小 , 而 焊 颖 处 材料 差异 和 轨道 表面 不 平 度 更 
容易 引起 天 线 轮 轨 间 作用 力 变化 , 该 应 力 的 突变 伴随 着 应 力 集中 , 对 天 线 轨 道 的 高 周 疲劳 有 
着 不 恨 影响 。 


本 文通 过 建立 的 QTT 轮 轨 模型 ， 考 虑 轨道 焊 颖 位 置 对 轮 轨 瞬 态 接触 的 影响 ， 研 究 了 接 
触 参数 和 轨道 面 形 貌 对 轮 轨 间作 用 力 的 影响 ， 主 要 分 析 结 果 为 : 〈1 ) 采用 层级 细 分 技术 对 模 
型 进行 前 处 理 ， 可 以 大 大 减 小 网 格 数量 和 计算 成 本 ， 且 模型 具有 较 高 求解 精度 。(2) 考虑 焊 
颖 材料 差异 时 ， 静 态 工 况 下 轮 轨 接触 总 变形 量 为 0.147 mm， 与 瞬 态 相 比 接触 变形 量 相差 较 
小 , 但 滚轮 经 过 焊 颖 时 轮 轨 间 Mises 应 力 发 生 剧 烈 变化 , 对 轮 轨 长 期 有 效 服役 带 来 不 利 影响 。 
G) 不 同 摩擦 系数 对 轮 轨 间作 用 力 影响 较 小 ， 主 要 影响 轮 轨 间 摩擦 应 力 ， 且 在 滚轮 经 过 焊 
颖 时 摩擦 应 力 增 大 约 71%。 对 大 载荷 天 线 轮 轨 间 作用 力 影 响 较 大 的 是 轮 轨 接触 面 形 貌 的 改 
变 ， 由 于 轮 轨 间作 用 力 大 载荷 、 小 集中 的 特点 ， 轨道 表 面 微小 不 平 度 特征 会 带 来 应 力 较 大 变 
化 ， 应 保证 天 线 轨道 表面 清洁 ， 时 常 检测 并 及 时 处 理 轨 道 表面 不 平 度 。 


Finite Element Analysis of Wheel-on-track Rolling Contact for 


Large Aperture Radio Telescope 


Shi Hongqiang!?, Xu Qian??^, Wang Na??^, Wang Zhaojun! 
( 1. School of Physics and Technology, Xinjiang University, Urumqi 830046; 
2. Xinjiang Astronomical Observatory, Chinese Academy of Sciences, Urumqi 830011; 
3. Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Urumqi 830011; 
4. Key Laboratory of Xinjiang Radio Astrophysics, Urumqi 830011 ) 


Abstract: Based on wheel-on-track rolling contact model which established by explicit finite 


element method, the stress variation about the wheel of 110-meter aperture radio telescope passing 
through the welding joint is analyzed. When modeling, not only the material difference at the 
welding joint, but also the influence of the flatness and friction coefficient at the track welding 
joint were considered. According to the model, when the wheel passing through the welding joint, 
the transient response of the force between wheel and track is analyzed. The hierarchical 
subdivision technology is used to preprocess the wheel-on-track model, and the results of static 
analysis show that the model has high precision. The transient analysis results show that, the 
Mises stress was increased by about 20% due to the track flatness, which may bring some hidden 
danger to the track's long-term reliable service. Therefore, it is necessary to inspection and 
maintenance the status of the track regularly to ensure the long-term reliable service 
characteristics of the track. 

Key words: Radio telescope; Azimuth track; Rolling contact; Transient response; Finite element 


method 
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